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1. Introduction and the idea of the method 

The mixed quark-gluon condensate, g (t^a^F^if)}, which by its definition measures the 
average interaction of the color- magnetic moment of a quark with the vacuum fields, plays 
an important role in QCD sum rules O (for a review see ||). Lattice measurements of 
this quantity in the chiral limit have been performed Q . Surprisingly, an analytic formula 
for this condensate has been derived only quite recently within the stochastic vacuum 
model of QCD g, ^| (for a recent review see [[?]]). In this paper, we will use the method 
of ref. ||] to evaluate the mixed condensate in the heavy-quark limit (i.e. for b, c, and t 
quarks), and for an arbitrary number of space-time dimensions. The method is based on 
the so-called world-line formalism ll, which is by now a well developed tool for perturbative 



calculations in field theory (for a review see [ 10 1 ) . Besides that, various attempts exist to 



apply the method to nonperturbative problems as well [11, H. Here, we will thus continue 
this line of research. 

The idea, which makes the world-line calculations possible analytically, is to transform 
the world-line integral to the one in an effective constant Abelian field, to be averaged over. 
The weight of the average, being prescribed by the form of the heavy-quark Wilson loop 
in the original confining theory, inherits therefore the information about that theory. The 
area of the minimal surface bounded by the quark trajectory C, S m i n (C), is < m~ 2 , where 
m is the current quark mass. For heavy quarks, S m i n (C) becomes smaller than the squared 
vacuum correlation length, Tr, since T g ~ lGeV" 1 0, [1|. (This length is defined as a 
distance at which the two-point gauge-invariant correlation function of F^ u 's falls off at 
least exponentially.) For such small-sized Wilson loops, one can write 



Smin — rX^liv (1-1) 
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Here, T,^ u = <f c x^dx v is the so-called tensor area, and "~" means "for nearly flat contours" , 
since eq. ( |1.1| ) is exact only when C is flat and lies in the (/i^)-plane. 

In the known confining theories, the following formulae for a small-sized Wilson loop 
exist 1 : 

• Stochastic vacuum model || 2 : 

(W(C))~iVexp(-CE^), (1.2) 

wh ^=803l)> <F 2 > ee <(F^(0)) 2 >; 

• Abelian-type theories with confinement: 

(W(C)) ~ N exp ( -a\l ^X 2 I . (1.:?) 




To the latter theories belong the 3d weakly coupled Georgi-Glashow model [14] and the 4d 
dual Abelian Higgs model. In this paper, we will consider a SU(iV)-generalization of the 
latter model, which was proposed for N = 3 in ref. and generalized to arbitrary TV in 



ref. lq| . This is just a dual superconductor model of confinement, based on the monopole 
condensation. With the logarithmic accuracy, the string tension of a short string in this 
theory reads a = 2tt(N — l)r] 2 ln/t, where rj is the vacuum expectation value of the dual 
Higgs field, and k is the Landau-Ginzburg parameter; ln/t 3> 1 in the London limit under 
study. 

Note that, in QCD at small distances, it is unlikely to have such a form of the Wil- 
son loop. Indeed, under the assumption that the Feynman-Kac formula, which relates 
the potential of a heavy qq-pair to the Wilson loop, can be extrapolated down to the 
distances smaller than T g , this would correspond to the linear next-to-Coulombic term in 
the potential. Such a term is however ruled out in QCD, as indicated by the pNRQCD 



calculations [17], the world-line calculations of the associated tachyonic gluon mass [18], as 



well as the lattice data [jl{|. As a support of this conclusion, it has been shown in ref. §§ 



that the heavy-quark condensate (?/>?/>) with the Wilson loop (1.3) is logarithmically UV 



divergent in 4d. Rather, the validity of the approximation (IT) for a heavy quark has been 
proven by demonstrating that, for the Wilson loop given by eq. (|1.2|) , the standard QCD 
sum-rules result for the heavy-quark condensate [^] is reproduced correctly. Neverthe- 
less, it is instructive to calculate the heavy-quark condensates, (4>ip) and g (ipcr^F^ip'j, 
in the case (1.3) as well. Firstly, because in the dual Abelian Higgs model [and its SU(iV)- 



generalization] a physical UV cutoff exists, that is the mass of the dual Higgs field, Mjj. 
Secondly, this calculation is interesting in order to find for the case (|1.3j), in the same 



way as for the case (1.2), the corresponding critical dimension of space-time. This can be 



x As in g §, (W(C)), which we consider in this paper, is a purely nonperturbative part of the full 
Wilson loop, i.e. the part, which can generate the terms starting from the linear one in the quark-antiquark 
potential. The Coulombic term corresponds to the one-gluon exchange, which yields the multiplicative 
renormalization of the loop. The corresponding renormalization factor is referred to the definition of the 



path-integral measure in eq. (2.1) below. 

2 This expression will be proven at the beginning of subsection 2.1. 
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defined as such a dimension d c , that, at d < d c the corresponding condensate is UV finite, 
it diverges as In ~ at d = d c , and as A d ~ dc at d > d c . Here is a table, which summarizes 



the values of d c : those for (ifiip) were found in 
found in the present paper. 



whereas those for g Njcr^F^ij)) will be 





{W>) 


g (tpcr^F^tp) 


case (1.2) 


6 


8 


case (1.3) 


4 


2 



It is further interesting to go beyond approximations ( |1.2| ) and ( |1.3| ) by accounting for 
the variation of the gauge field, which is significant at the distances of the order of T g . In 
this paper, we manage to do so and find, within the stochastic vacuum model, an explicit 



correction to the formula (1.2), in case when C is a circle. This allows us to estimate (at 
least parametrically) corrections to (V'V') an d g {^cr^F^if^ , which are produced by the 
variation of the field. In particular, as far as their signs are concerned, we find that both 
corrections are negative, i.e. they diminish the absolute values of the condensates. 

The paper is organized as follows. In the next section, the mixed condensate will be 



derived in both cases, (|1.2| ) and (1.3). In section 3, these results will be generalized to 
finite temperatures. In section 4, corrections emerging due to the variation of the gauge 
field will be evaluated. In section 5, a brief summary of the obtained results will be given. 



2. Mixed quark-gluon condensate at zero temperature 



To evaluate the condensate g {^cr^F^yijj^ , we will use the world-line representation for 
this quantity § and write explicitly the integral over the Grassmann functions ip^s, which 
describe spin degrees of freedom (cf. refs. || |i~0f): 



v 



o 



— — I m — 7 M 



dTe~ m2T x 



A" 2 



Vx,. 



£>W exp 



dr 



x^(0)=a; M (T) 



</v(o)=-Vv(T) 



x J dTa^a Xp ti(F f , v (x(0))F Xp (x(T))){Pexp 



T 



dr (A^x^ - ipfj.ipuFf.u) 



(2.1) 



Here, V is the four-volume occupied by the system, A stands for the UV momentum cutoff, 
= mTfij lv] with 7^'s being the Euclidean gamma- matrices, the average (...) is defined 
with respect to the gluodynamics action in the Euclidean space-time, j J d A x (i 7 ^) 2 , where 
a = I, ... ,N 2 — I and F* v = d^A^ - d u A a ^ + gf abc A b ^Al stands for the Yang-Mills field- 
strength tensor. Next, A^ = A^T a with T a 's standing for the generators of the SU(iV)- 
group in the fundamental representation, [T a ,T*] = if abc T c , tr T a T b = \5 ab . Finally, on 



-3- 



the r.h.s. of eq. (|2.l| ) 4 = 2(2s + 1), where s = 1/2 is the spin of a quark, and the factor 
(2s + 1) counts the number of spin degrees of freedom 10|. Rather, an extra factor 
2 appears from the expression for the mixed condensate as a variational derivative with 
respect to the coupling g (which is made x-dependent for a while Q) of the averaged quark 
propagator, which is defined at a closed path: 



g (tpo-^Ffj,^) = tr (g(x)ff^ v F^(x)S(x,x\A^)) = 2tTj-^-(S(x,x\A^)), where 

S (x, y\Al) = (m + 7^)^, D„ = d„ - igA*T a . 

As in ref. ||, by virtue of the operator of the area derivative, 5/5s PjL ,(x(t)), all the 
gauge-field dependence of this expression can be reduced to that of the Wilson loop 
(W{C)) = (trVexp (ig J Q T drA^j J). Indeed, 



tr(^(i(0))^(i(T))> (Pexp 



ig I dr (A^Xfj, - ip^F^) 
o 



g 2 5s pu (x(0))5sx p (x(t)) 6Xp V " J "" Ss pu (x(t)) 



T 



drip u A 



(W(C)) 



Further, by splitting in the standard way the coordinate x p (t) into the position of 
a trajectory and the relative coordinate, one gets from the empty integration over the 



position the factor V, which cancels with 1/V on the r.h.s. of eq. (2.1). The relative 



coordinate will be denoted below as z^{t). 



2.1 Stochastic vacuum model 



To start with, let us present a derivation of eq. ( p..2j ) from the stochastic vacuum model. 
This model suggests the following parametrization for the nonperturbative part of the 
two-point correlation function of gluonic field strengths: 



1 



NxN 

N 



■M< (S^8 vp - 8„ p 5 vX )D ((x - x') 2 ) + 



+^ K {{* ~ x')x5 up - (x - x') p 8 vX ) + d x v {{x - x') p 8^ -(x- x') x 5 pp )] D l {{x - x'f) \. 



Here, <fr xx i = i V exp (ig f®, dz p A p (zj) is a phase factor along the straight line, which goes 
through x' and x, and 

<F 2 > 

N = 4n[D(0) + D 1 (0)] (2 ' 2) 

is the normalization constant. The corresponding expression for the Wilson loop reads (see 
e.g. pi] for the details): 



(W(C)) = exp< 



8iV 



2 / dff^lx) [ da^ix^D ((x - x') 2 ) + 
Jt,(c) Jh(c) 
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+ 



j> d XfJL j> dx'^G ((x - x') 2 ) 



G(x 2 



dtDxit). 



(2.3) 



Let us now perform a Taylor expansion of this expression at the distances |x| < T g of 
interest. In the leading approximation, we have 



G(x 2 ) ~ G(0) -x 2 D 1 (0), D(x 2 )~D(0), 



(2.4) 



that yields 



j d Xfl j> dx'^G ((x - x' f) ~ 2£>i(0)Ej 1/ , 
da^(x) / da„ v {x')D ((x - x') 2 ) ~ D(0)E 2 . 

2(C) JE(C) 



Taking into account the value of the normalization constant A/ - , eq. (|2.2| ), we arrive in this 
approximation at eq. ( |l.2|) . In section 4, by finding the next term of the Taylor expansion, 



we will obtain corrections to eq. (1.2) and to the quark condensates produced by the 
variation of the gauge field. 



Further, the square of the tensor area in the exponent in eq. (1.2) can be linearized 
by using the Hubbard-Stratonovich transformation with a constant-valued antisymmetric- 
tensor field B, 



iVexp(-CE 



(IV ) 



N / TT f + °° _<L 



exp 



■ 2 Biw^j/w j - 



is a half of the number of off-diagonal components of B^ u in d dimen- 
sions. The operator 5/5s^ v can now readily be applied, since, due to the Stokes' theorem, 
B^u^^u = /we) da^yB^y. In particular, since the field B^ v is space-time independent, 
the operator 8/8x^(0) in eq. (2.1) yields zero. The world-line integral becomes the one 



where n 



in a constant Abelian field, and the result is nothing but the Euler-Heisenberg-Schwinger 
Lagrangian (see e.g. pi 



z M (0)=« M (T) 



</v(0)=-</y(T) 



x exp 



1 1 ■ i 



1 



[T 2 abcot(aT] coth(6T) - l] . 



(4irT) d / 2 



(4vrT) d / 2 

Here, we have in the standard way subtracted the part of the integral without the field, 
in order to respect the obvious normalization condition g (^(J^yF^yil) 



0. We have 



also used the standard notations (lCj 
o 1 



E 2 - H 2 + a/(E 2 - H 2 ) 2 + 4(E • H) 2 



- 5 - 



E 2 - H 2 ) + y (E 2 - H 2 ) 2 + 4(E • H) 2 



with E = i (B41, -642,-643) and H = (B23, — -B13, B12) being the electric and magnetic fields, 
which correspond to the field-strength tensor B^ [By = e^H^i B^ = —iEi). 

Next, due to the factor e~ m T , this expression can be expanded in powers of T: 

T 2 ab cot (or) coth(6T) - 1 = 
= y (b 2 ~ a 2 ) + O (T\E ■ H) 2 ) = ^Y, B ^ + ° ( T4 (^) 2 ) ■ ( 2 - 5 ) 

a<0 

Note that, in terms of the fermionic determinant, this approximation corresponds to re- 
taining only the diagram with two external lines of the field B^ v in the expansion of the 
logarithm. The parameter of the expansion can actually be estimated accurately. To do 
the estimate, let us set d ~ 1, having in mind the values of d c given by the table in sec- 
tion 1. Then, since T ~ mT 2 and B 2 U ~ g 2 (-F 12 ) /N, the parameter of the expansion is 

jj m 4 I , i.e. the heavy-quark limit is implied in the sense of the inequality 

^? < 1. (2.6, 

For N = 3, the expansion ( p.5| ) breaks down for u, d and s quarks, but holds for c, b and t 
quarks under study. Note also that the parameter of the expansion is O [N ) at large N. 

As follows from eq. ( |2.5[ ), the adopted approximation reduces the range of the B^ v - 
integration to the interval (— 1/T, 1/T) ~ (—m 2 ,m 2 ). However, the characteristic values 



of yj B 2 U , being of the order of O y\f g 2 (F 2 ) /NJ , are smaller than m 2 , as long as the 
expansion Q2.5D converges. This fact enables one to retain the infinite range of the B^ u - 
integration with the exponential accuracy. Apparently, this is true only for the case ( |l.2| ) 
under study, when the measure of the -B^-integration falls off so rapidly. Rather, in 
case (|l.3[), which will be considered in the next subsection, the measure will be shown to 



fall off only as some power of B 2 U , and the infinite range of integration over B^ u cannot 
be retained anymore. 

The appearing integral over B^ can then be readily calculated: 

(2.7) 

Bringing all the factors together, taking the limit A — ► 00, and noticing that <7 2 „ = ^, we 
have 

,- 128n(f + l)iV a [°°dTe- m2T 
9 (^^F^) = \ 2 ' mC 2 / - d 




3(47r) rf / 2 J T l-3 ' v ' ; 

For the case d < 8, where the integral is convergent, we obtain the final result (valid at 
d > 2): 
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where 'T" henceforth stands for the Gamma-function. As well as the standard heavy-quark 
condensate {ipip), this expression scales at large N as 0(N). At d < 6, when (^V) is a l so 
UV finite H , we obtain for the ratio of the two condensates 3 



2_ 9 ff(jV^) (I + 1) (3 - |) ff 2 (F 2 ) 



In particular, 



[ijjTp) nN 
2 9 2 (F 2 ) 



in,-. 



l 0| d =4 3 ^2 • 

2.2 Dual superconductor model 

In this case, one needs an additional integration over an "einbein" parameter, in order to 
get rid of the square root in the exponent in eq. (1.3): 

{W (C)) = N r*Lexp(-\- 



\ 8A 



o 



Using further again the Hubbard-Stratonovich transformation and integrating over A, one 
gets the following expression for the Wilson loop §§: 

The weight of the ^^-integration in this case is therefore only polynomial and not Gaus- 
sian. Let us adopt again the approximation (2.5), where the expansion parameter now 
becomes a/m 2 , that establishes an upper bound on the value of a. Taking for instance the 
minimal value of m, m c = 1.1 GeV, and parametrizing a = Acioo, we obtain A < 6.7. 
The .B^-integration then takes the form [cf. eq. (j2.7|)]: 

B uv B\ n 2 

rL^ = 




1 (,2 + L ) JO (l + X 2 

Denoting t = m 2 T, we arrive at the following intermediate result: 

2 3-d NT (n±|) a A m d-7 J ^ _ t °i dxx n+3 



9 {^^F^lp) = 2 d+1 / dtt 



2C 



n+1 



K2 > (m/A) 2 v ^ 

Clearly, the integral over x here receives the dominant contribution from the region of t 
around (m/A) 2 . For such t's, the x-integral approximately equals 3(7^-) j and we arrive 



3 As in ref. ffl], this ratio is defined with the factor 2, in order to bring it in accordance with an analogous 
quantity used in QCD sum rules. The reason is due to the different definitions of <r M „: <7 M „ = s* [7^,71/] in 
QCD sum rules, while cr pi/ = jj[7h,7k] in [§] and in the present paper. 
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at the following result: at d = 2, the condensate diverges logarithmically in the UV region, 
namely 

/ - \i AN am A 

ff(^F^)| d=2 = — (3^r ln -> ( 2 - n ) 

while, at <i > 2 it diverges as A rf ~ 2 , i.e. polynomially: 

2^r(«+l)iV ( 7mA' i - 2 

= - 9(i i^ r(s) ■ < 2 - 12 > 

Setting for A the mass of the dual Higgs boson, Mh, we can obtain a lower bound for 
the value of the Landau-Ginzburg parameter k = Mn/My, where My is the mass of the 
dual vector boson. Indeed, the condition A>m yields 

Tfl 

ln/t»ln . (2.13) 

My y 1 

On the other hand, it can be shown [^jj that My 1 is the distance, at which the two-point 
correlation function of F^s in the dual Abelian Higgs model, or its SU(A r )-generalization, 
falls off. Therefore, My should be of the order of T~ . For c and b quarks, this value of 
My leads to the condition hiK 3> 1, that is the standard condition of the London limit in 
the type-II (dual) superconductor [p4f| . Rather, for t quark, whose mass is about 170 GeV, 
this condition should be replaced by Inn 3> 5.1. Therefore, it makes sense to speak about 
the heavy-quark condensates in this model, provided the condition of the London limit is 
fulfilled in a certain, m-dependent, way. 



3. Finite-temperature generalizations 



In this section, we will generalize the above-obtained 4d-results to the case of finite tem- 
peratures. In case (|1.2|) , we will adopt the QCD-value for the temperature of dimensional 
reduction, T<j. r . ^ 2T C (see e.g. J?]]). At temperatures smaller than Td. r ., it is important to 
account for the antiperiodic boundary conditions for quarks, that can be done upon the mul- 



tiplication of the zero-temperature heat kernel by the factor [25] 



1 + 2 E("l) n e 
n=l 

Here, (3 is the inverse temperature, and n is the number of a Matsubara mode. Inserting 
this factor into eq. ( |2.8D and calculating the integral over the proper time, we have for the 
case flL2|) : 

rf) (T) = G(T) = 



= G(0) 

where G(0) = 
(3m » 1. At T 



1 + {(3m) 2 Y2(-l) n n 2 K 2 (Pmn) 



n=i 



G(0) [l - 072(/3m) 3 / 2 e-^ 



i 



44t^ — Wm? — • This result is valid at T < T c , where, for heavy quarks, 
= T c , the chromoelectric condensate evaporates, and one should substitute 



2g 2 ({B^) 2 ) instead of g 2 (F 2 ). Note that the lattice data @ also confirm that, at T = T c , 
g 2 {F 2 ^ drops by a factor of the order of 2. With this substitution, the obtained result is 
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valid up to Td r , since (3m > 1 up to this temperature. At T > T^_ T _, the sum over Mat- 
subara frequences disappears, and one is left with a 3d theory with the choromomagnetic 



condensate. The mixed condensate in this theory is given by eq. (2.9) at d = 3. As follows 
from the elementary dimensional analysis, which relates fields and couplings in this theory 
to those in the original 4d one, the mixed condensate in the original QCD at T > Td. r . 
stems from this result upon the multiplication by T. It reads 



G(T) 



5T {g 2 ((B-) 



96vriV 



Let us now consider the dual Abelian-Higgs-type theory at finite temperature. Since 
M H (T) oc y/a(T), where a(T) ~ (l - j)" at T -> T c - 0, we have 



M H {T) M H (0) 



m 



T 

1 - — 



v/2 



> 1. 



Since Mff(0) = K,My(0), we obtain a new condition for k [cf. eq. (]2 . 1 3|) ] : 



In k 3> In 



rn 



My(0) 



In 1 



T 



(3.1) 



Here, v should be set 1/3 and 0.63 |£7J for the 1-st and the 2-nd order phase transition, 
respectively. The obtained formula defines the law, according to which the constraint 
imposed on k should become more severe when T approaches T c . Provided this constraint 
is obeyed, we then have at T < T c : 



G(T) = G(0) 
Here, G(0) = 



1 + 



8mT 



E 



-l) n 



n 



K\ (f3mn) 



G(0) 



1 



M%(T) 



5Nm ^y>^' T , and we have assumed that, as well as in QCD, T c -C m. 



4. Accounting for the variation of the gauge field within the stochastic 
vacuum model 



In this section, we will estimate corrections to (V'V') an d tuq at T = 0, which can be obtained 
within the stochastic vacuum model by treating the gauge fields as varying, in the leading 
approximation. To this end, we need to find the next term in the Taylor expansion of 
eq. (|2,3j). This is only possible if we specify the form of the functions D and D\ at the 
distances \x\ < T g of interest. In ref. [28|, it has been argued that, at such distances, 



these functions should have a Gaussian shape. This follows for instance from the formula 



dD(x 2 ) 
dx 2 



x=0 



oc gf 



abc 



fpa 77b 77c 

/us v\ \fX 



which, due to the finiteness of the triple condensate, 

means that D should be an analytic function of x 2 at short distances. We therefore take 
the functions D and D x in the form D(x 2 ) = D(0)e- x2/T s, Di(x 2 ) = Z?i(0)e~ x ' 2/T 9 . The 

D(0) 



extra terms, one then gets in eq. (|2.4j ), read AG(x 



£>i(0) 



'1 

(x 2 ) 2 , AD{x 2 



9 



-x 



To 



estimate the effect, which these corrections produce on the condensates, we will assume for 
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simplicity that C is a circle of the radius R. We then obtain very similar expressions for 
the corresponding integrals: 



jf d Xfl jf dx'^AG [(x - x'f) = -IG^^i? 6 , 
2 f da^(x) [ da^{x')AD ((x - x'f) = -16vr 2 ^i? 6 , 

JT,(C) JS(C) 1 g 

where in the last case we have used the formula da /ll/ (x) = ^{x^dx v — x v dx^). Combining 
together the leading and the next-order results, we obtain 

9 2 (F 2 ) 



(W(C)) ~ iVexp -CZ* + a{T,l u f 2 



2 3 / 2 7riV(d 2 -d)T}.' 



fi (s 2 ) 3/2 

where we have used that R = 2 ^J 2 tt z ' ra tio of the correction to the absolute value 

ct(£ 2 ) 3 / 2 

of the leading term is — -^r^ — ~ ■ This means that the correction is anyway small 

as long as the size of C, that is the quark's wavelength 1/m, is smaller than T g — the 
approximation we work in. To perform the calculations, we will, however, need the validity 
of a more fine inequality, which will be derived in a moment. 

Apparently, the correction can be analysed if we get rid of the power 3/2. This can 
be done by noticing that e Q ( s ^ 3/ ' is nothing but the saddle-point value of the integral 
c Jq 00 dX exp (— 2^3- A 3 + ^fiu) ■ The normalization constant c here is fixed by the condition 

1 = c J d\e 27q 2 , which yields c = p(i/3) ' h na l expression for the Wilson loop 
therefore reads 

According to this formula, the requirement for the next-order result to be small with respect 
to the leading one implies that the characteristic values of A are to be much smaller than 
C. Parametrically, these characteristic values are of the order of a 2 / 3 , that leads to the 
condition N ^ S> 1. This inequality apparently breaks down at large d, that does not 
matter anyway, since the condensates diverge at d > d c . At d < d c , we therefore have, 



together with the heavy-quark limit condition (2.6), the following sequence of inequalities: 



1 / N \ l ' A 



m~{g 2 {F 2 )) <<T9 ' 

It clearly tells us that the quark wavelength is smaller than that of the gauge field, which 
itself is much smaller than the vacuum correlation length. The first inequality means that, 
inside C, the gauge field is effectively constant, while the second one means that the field 
is varying significantly at the scale of the vacuum correlation length. 

The last step in evaluating the desired corrections to the condensates essentially 
amounts to substitute C — > C — A, C 2 — > C 2 — 2CA and integrate over A. We obtain 

(#} new = (V^> (1 - 0» 9 (^^^F^) new = g (^0>iV0) C 1 - 2£)> "ig new = mg(l - £), 

(4.1) 
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where 



(2/q) 2 /3 ( 4 x3 ^ _ 3 • 2 4 / 3 • T(2/3) \ N{d 2 - dy 1/3 

'i 1 r- : r} 



r(i/3)c 7 dAAexp V 27« 2 ' v ;~ 7r 2 /3.r(i/3) 



(4.2) 



and (VV') = 3(4^-)d/2 — 9 2 i^ 2 ) * s the value of the heavy-quark condensate at 2 < d < 

6 H without the correction. This result gives an idea of how the variation of the gauge 
field affects the condensates in the leading approximation. In particular, we see that the 
corrections diminish the absolute values of both condensates, as well as the value of to 2 ,. 

5. Conclusions 

We have calculated the mixed quark-gluon condensate, in the heavy-quark limit, at zero 
and finite temperatures. This has been done in the models of stochastic vacuum and 
dual superconductor, where the nonperturbative parts of the Wilson loop at the distances 



smaller than the vacuum correlation length are given by eqs. ([h^) and (1.3), respectively. 
In the stochastic vacuum model, the condensate is UV divergent at d > 8, while in the dual 
superconductor it diverges right at d > 2. This divergence can be regularized by attributing 
to the UV cutoff a physical meaning of the inverse thickness of a "short string", that is 
nothing but the mass of the dual Higgs field. We have also found the constraints, which 
the Landau-Ginzburg parameter should obey at zero and finite temperatures, in order for 
such an interpretation to be possible [cf. eqs. (|2.13|) and (|3.1[) 1- At zero temperature, the 
obtained expressions for the condensate are given by eqs. (§J), (pll| ), and (|1|), while 
the finite-temperature generalizations are discussed in section 3. 

We have further evaluated corrections to (dip) and g {^o^F^t^ , which appear within 
the stochastic vacuum model, due to the variation of the gauge field at the scale of the 
vacuum correlation length. To this end, we have taken into account that, at the distances 
smaller than T g , the two-point correlation function of gluonic field strengths has a Gaussian 
shape. The calculation has been done in the most simple case, namely when the quark 
trajectory is a circle. Despite this simplification, the obtained formulae, eqs. (|4.1| ) and (|4.2j), 
yield the parametric form of the corrections and show that these diminish the absolute 
values of the condensates, as well as the ratio of the mixed condensate to the standard one. 

The world-line approach to the calculation of heavy-quark condensates, proposed in 
ref. H and in the present paper, has several positive features. First of all, as we have 
seen, it allows to calculate the condensates in various models of the QCD vacuum, where 
a certain nonperturbative (^-potential at small distances is generated. In particular, the 
stochastic vacuum model, where this potential is quadratic, reproduces the 4d-value of 
{'ipip), which stems from QCD sum rules. We have also seen that the proposed approach 
naturally allows for the finite-temperature generalizations. Since the influence of thermal 
effects on the condensates is an important problem, it would be very interesting to have 
the lattice data on m^T) not only in the chiral limit p9|, but also in the heavy-quark one, 
and to compare them with our results. We have further seen that the world-line approach 
is appropriate for a systematic analysis of corrections to the condensates, which originate 
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from the variability of the gauge field inside the quark trajectory. Last but not least, 
our approach can also be used for the calculation of the higher mixed heavy-quark-gluon 
condensates of the form g n ^(a^ u F^ u ) n ip^, where n > 2 is an integer. 

Note in conclusion that recently |3(|, i^if) has been obtained as a function of m/T, by 
doing calculations on the gravity side of the AdS/CFT correspondence. The large- m part 
of that result qualitatively agrees with the 1/m fall-off, which is predicted by the stochastic 
vacuum model and QCD sum rules and disagrees with the behavior mln^, which stems 
from the dual superconductor model ||. It would be very interesting to have analogous 
results for g ( ^a^F^ip) as well, and to compare those with the predictions of the present 
paper. 
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